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We study the properties of the ++ glueball at finite temperature using SU(3) 
quenched lattice QCD. We find a significant thermal effects near T c . We perform 
the x 2 fit analyses adopting two Ansatze for the spectral function, i.e., the conven- 
I tional narrow-peak Ansatz and an advanced Breit-Wigner Ansatz. The latter is an 

O ■ extension of the former, taking account of the appearance of the thermal width at 

T > 0. We also perform the MEM analysis. These analyses indicate that the ther- 
mal effect on the glueball is a significant thermal- width broadening r(T c ) ~ 300 
MeV together with a modest reduction in the peak center Auiq(T c ) ~ 100 MeV. 



Near the critical temperature T c , the QCD vacuum begins to change its 
structure. For instance, the string tension reduces and the spontaneously 
broken chiral symmetry begins to restore partially. Since hadrons are com- 
posites of quarks and gluons, whose interaction depends on the proper- 
ties of the QCD vacuum, it is natural to expect the structural changes of 
various hadrons. In this direction, there are a number of studies based 
on QCD motivated effective models, which predict the mass reductions of 
charmonium and light qq mesons near T c . 1,2 In BNL-RHIC experiments, 
these changes are considered as important precritical phenomena of QCD 
phase transition. In addition, lattice QCD Monte Carlo calculations have 
been performed to examine the pole-mass shifts of hadrons at the quenched 
level. 3,4 The ++ glueball is one of such hadrons, which are expected to ex- 



*The numerical calculation has been done on NEC SX-5 at Osaka Univeristy. 
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hibit structural changes near T c . Studies based on the effective models 5,6 
suggest a large mass reduction of the glueball as a natural consequence of 
the sizable difference between the ++ glueball mass ttiq — 1500 — 1700 
McV and the critical temperature T c = 260 — 280 MeV. In this paper, we 
investigate the glueball properties at finite temperature by means of lattice 
QCD Monte Carlo calculation. 7 ' 8 ' 9 

We consider the temporal correlator G(t) = (4>(t)4>(0)) of the glueball 
operator 4>{t). The simplest construction of the glueball operator is the pla- 
quette operator, which, however, is known to have only a negligible overlap 
to the lowest-lying ++ glueball. To enhance the overlap, the smearing 
method provides us with a systematic procedure by extending the spatial 
size of the operator. In this paper, we adopt the APE smearing. 8 ' 10 The 
continuum expression of the smeared glueball operator is obtained in the 
Coulomb gauge as 

J ($^3 -p ( J1 1 ) <w<w, (i) 

where Gfj is the field strength tensor, and p <E R controls the size of the 
spatial extension. (For detail, see Refs. 8.) By choosing p appropriately, we 
can maximize the contribution from the lowest-lying glueball in the glueball 
correlator G(r). 

We next consider the spectral representation of G(t) = (0(r)</>(O)) as 

/>oo 

G(r) = / dwK(T,w)A(u), (2) 
Jo 

with (3 = 1/T, K{t,uj) = cos g ^h(£/2) T)) ' and denotes the spectral 

function with its spatial momentum projected to zero. At low temperature, 
peak positions of A(ui) provide us with the pole masses of hadrons. Hence 
the simplest parameterization of A{ui) for the x 2 fit analysis is given as 7 

A{u) =C{5(u-m) -6(u + m)}, (3) 

where G and m are used as the fit parameters, corresponding to the strength 
and the pole mass, respectively. To adopt this Ansatz, the peak should be 
narrow enough. We will refer to Eq. (3) as the " narrow-peak Ansatf . 
At T > 0, however, through the interaction with the thermally excited 
particles, the thermal width is generated even for stable particles. As T 
gets larger, the effects of the thermal width become more important. To 
respect its existence, we adopt the following parameterization as 8 



A(u) = C {S r (cj - w ) - 5 r (uj + w )} , 



(4) 



February 1, 2008 4:25 WSPC/Trim Size: 9in x 6in for Proceedings 



ishii 



3 



where 5r(u)) 




. C, u>o, and T are used as the fit parameters 



corresponding to the strength, the peak center, and the thermal width, re- 
spectively. We will refer to Eq. (4) as the u Breit- Wigner Ansatf . Note that 

lim dr(u)) = S(u)). Hence, the Breit- Wigner Ansatz is a natural extension 
r^o 

of the narrow-peak Ansatz. 

For these \ 2 fit analyses, proper Ansatze for A(u>) have to be provided. 
However, in the very vicinity of T c , and above T c , there may appear more 
complicated structures. In this sense, it is desirable to obtain A(u>) directly 
from G(t), which, however, is known to be an ill-posed problem numerically. 
Recently developed Maximum Entropy Method (MEM) 11,12 can deal with 
this inverse problem adopting the following Shannon-Jaynes entropy as 



where m{uu) is a real and positive function referred to as the default model 
function. m(oj) is required to reproduce the asymptotic behavior of A(oj) 
as w — > oo. 

We use the SU(3) anisotropic lattice plaquette action 13 as 



S G = ReTr{l-P y ( S )} + ^^ReTr{l-P 44 ( S )}, 



where -P^(s) G SU(3) denotes the plaquette operator in the \irV- plane. 
The lattice parameter and the bare anisotropic parameter are fixed as 
/3i at = 2N c /g 2 = 6.25 and 7g = 3.2552, respectively, so as to reproduce the 
renormalized anisotropy as a s /at — 4. The scale unit is introduced from 
the on-axis data of the static inter-quark potential with the string tension 
\fa = 440 MeV. The resulting lattice spacings are given as a^ 1 = 9.365(66) 
GeV and a^ 1 = 2.341(16) GcV. The critical temperature is estimated as 
T c ~ 280 MeV from the susceptibility of the Polyakov loop. 8 We generate 
5,500-9,900 gauge configurations to construct the glueball correlators, 7 ' 8 
where statistical data are divided into bins of the size 100 to reduce pos- 
sible auto-correlations near T c . To enhance the low-energy spectrum, we 
adopt an appropriate smearing corresponding to p ~ 0.4 fm in Eq. (1). 

We present the numerical results. Whereas the narrow-peak analysis 
suggests that the thermal effect on the glueball is the pole-mass reduction 7 
of about 300 MeV, the advanced Breit- Wigner analysis suggests that it is 
a thermal-width broadening 8 of about 300 MeV together with a modest 
reduction in the peak center of about 100 MeV. The numerical results of 



S = 





s,i<j<3 



s,i<3 



(6) 
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the x 2 fit analyses with the Breit-Wigner Ansatz are shown in Fig. 1. 
C0 [MeV] r [MeV] 
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Figure 1. The peak center u>o(T) and the thermal width T(T) of the ++ glueball 
against temperature T, obtained from the x 2 fit analysis with the Breit-Wigner Ansatz. 
The vertical dotted lines indicate the critical temperature T c ~ 280 MeV. 

We now turn to the result of MEM analysis. As the default model 
function, we adopted the perturbative expression of the smeared glueball 
correlator up to O(ag) as 9 

m(w) = Nuj 4 exp {-(up) 2 /A) , (7) 

where the normalization factor N is determined from 1 = J °° (koK(r = 
0,u>)m(u>) so as to mimic G(t = 0) = 1, which we adopt as the normaliza- 
tion of the correlator. Following the numerical procedure given in Ref. 12, 
we reconstruct A(lo). The results at various temperatures are shown in 
Fig. 2. We see that the peak tends to become broader with the increasing 
temperature, 9 which is consistent with the results of the Breit-Wigner fit 
analysis. 

To summarize, we have studied the properties of thermal glueballs be- 
low T c by using SU(3) anisotropic lattice QCD at the quenched level. We 
have performed \ 2 fit analyses adopting two Ansatze for the spectral func- 
tion A(u>), i.e., the conventional narrow-peak Ansatz (5-function type) and 
the advanced Breit-Wigner Ansatz so as to respect the appearance of the 
thermal width at T > 0. Whereas the former has suggested that a signif- 
icant pole mass reduction of the ++ glueball of about 300 MeV near T c , 
the latter has suggested that the thermal effects on the glueball is rather 
a significant thermal-width broadening of about 300 MeV together with a 
modest reduction in peak center of about 100 MeV. We have also performed 
MEM analysis. We have seen that the peak of the ++ glueball tends to 
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Figure 2. MEM results for the ++ glucball spectrum at various temperatures below 
T c . 

become broader with the increasing temperature below T c . This is consis- 
tent with the results of the advanced Breit-Wigner fit analysis, suggesting 
the importance of the thermal width. Experimentally, there are two glue- 
ball candidates, 14 i.e., /o(1500) and /o(1710), both of which have rather 
narrow decay width T ~ 100 MeV. Hence, the experimental observation 
of thermal-width broadening of 300 MeV would be possible in the future 
experiments at RHIC. For more realistic comparison with experiments, it 
is desirable to estimate the effects of the dynamical quarks on the thermal 
properties of the glueball. 
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